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I present preliminary results on the search for h c in its rj c j and J/ipiv decay modes. We observe an excess of 
■q c -y events near 3526 MeV that has a probability V ~ 0.001 to arise from background fluctations. The resonance 
parameters are M = 3525.8±0.2±0.2 MeV, F < 1 MeV, and 10.6±3.7±3.4(&r) < Y^ p B Vc ^ < 12.8±4.8±4.5(6r) eV. 
We find no event excess within the search region in the J/ipn° mode. 



1. INTRODUCTION 

Charmonium states have been succesfully stud- 
ied in pp annihilations, where all states can be 
formed, and detected in their decay to electro- 
magnetic final states. E760 and E835 measured 
precisely masses and widths of rj c , J/ip, ip{2S) 
and Xcj states as well as other properties for these 
states, e.g., T 77 and B p p 1 . 

The h c is so far the most elusive charmonium 
state below DD threshold |5] . It has not been ob- 
served yet in tf>{2S) decays, where it could be pro- 
duced by sequential radiative transitions through 
the Xc2 (El followed by Ml), or by I spin violating 
7T° transition, nor in B decays whose branching 
ratios to h c could be as large as O(10 -4 ) [3]. 

Observation of the ft. c (l 1 Pi(l H )) will allow to 
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determine the splitting between singlet P and the 
spin- weighted average mass of the triplet P states 
(Xc.o.g = 3525.30 ±0.07MeV gj). In one-gluon 
potential model, this splitting is zero. Corrections 
are small and the mass difference between the h c 
and the Xc.o.g is in general predicted to be at most 
a few MeV \Q , thus it is important to measure the 
h c mass to better than w 1 MeV. 

The h c is expected to be narrow (< lMeV in 
width), to have a dominant El transition to rj c j, 
and large branching ratio to light hadrons 6 . 

Our study of XcO has shown that El transitions 
of triplet P states are in excellent agreement 0j 
with the predicted scaling as the third power of 
photon momentum. Assuming that the h c has the 
same radial wave function as the other P states, 
this would imply T(h c — > i] c j) « 600 keV for an 
h c mass close to the Xc.o.g- 

In 1992, E760 (our former experiment) re- 
ported observation of a structure in the cross sec- 
tion pp^J/ijjw (an I- violating mode) close to the 
Xc.o.g. interpreted as the /i c [7|- Neglecting in- 
terference with continuum pp—>J/ipir°, resonance 
parameters were determined as M = 3526.2 ± 
0.15±0.2 MeV; T < 1.1 MeV (90% CL); and (1.8± 
0.4) • 10~ 7 < B{pp)B{J/ipiT ) < (2.5 ± 0.6) • 10~ 7 . 
The probability for such structure to arise from 
background was estimated to be 1/400. 

Observation of the h c is one of the principal 
objectives of experiment E835 at Fermilab. 
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2. EXPERIMENT E835 

Charmonium states are studied by a scan of 
the pp annihilation cross section for exclusive final 
states at different center of mass energies {Ecu)- 
An excess of events at any value of E^ M over the 
background measured on a broader Ecu range 
signals the formation of a resonance. Resonance 
parameters are then determined with precision up 
to 100 keV on masses and widths. The experi- 
ment is designed to observe charmonium states in 
their decays to electromagnetic final states such 
as J/ijiX — > e + e _ X and 7^(7) — > 77(7). 

E835 is a major upgrade of E760 and is de- 
cribed in detail in [8]. The detector is a non- 
magnetic, large acceptance, cylindrical spectrom- 
eter, covering the complete azimuth (0) and from 
2° to 70° in polar angle {&). It consists of a lead 
glass electromagnetic calorimeter divided into a 
barrel and a forward sections; the inner part of 
the barrel section is instrumented with a multi- 
cell threshold Cerenkov counter, for electron de- 
tection, three concentric scintillator hodoscopes, 
and a tracking system to measure charged parti- 
cles. To withstand the ~ 5 MHz continuous in- 
teraction rate, all channels are instrumented with 
multi-hit TDCs. 

The required high luminosity is achieved by 
a H2 jet target intersecting the p beam in the 
Fermilab Antiproton Accumulator. The beam is 
decelerated from the accumulation energy to the 
value appropriate for the formation of each reso- 
nance. The density of the target is continuously 
adjusted to compensate for beam loss keeping the 
instantaneous luminosity constant at ~ 2 x 10 31 
cm s . The stochastic cooling keeps beam mo- 
mentum constant with a Ap/pw 2 x 10~ 4 , com- 
pensating for energy losses in the target. 

The integrated luminosity Li nt for each energy 
setting is measured to < 2% by counting recoil 
protons from pp scattering at 9 w 90°. 



3. SEARCH FOR h c -> n c j and J/^jn 

E835 took data in 1996/1997 (E835-I) for an 
integrated luminosity Lm ~ 140 pb , and again 
in 2000 (E835-II) for L int « 110 pb" 1 . 



We search for h c in the following reactions 

pp — > h c —> J/ipTT° — ► e + e~77 (1) 

pp-*h c -+ i] c j -> 37 (2) 

The scan for the h c was based on assump- 
tions for the mass (close to Xc.o.g.), total width 
(<1 MeV) and expected yields at resonance 
peak: 3-8 ev/pb _1 (above a continuum yield of 
1=3 2 cv/pb -1 ) for reaction JIJ and < lev/pb _1 
for reaction UJ. 

We have taken data for about 215 pb -1 in the 
range 3300 < E C m < 4400 MeV. Approximately 
70 pb -1 were spent in a fine scan of the Xc.o.g. 
region 3525.2 < E CM < 3527.2 MeV and 20 pb^ 1 
a coarser scan of 3520 < Ecu < 3540 MeV (see 
Fig-EJi while the remaining data are used to mea- 
sure background. Data taken at Xd an d Xc2, on 
which we performed repeated scans which will 
provide new measurements of their masses and 
widths, provide also clean J/tpj events to monitor 
efficiencies and nearby background measurement 
for reaction (J2J) on cither side of the Xc.o.g- 

3.1. Search for pp — > h c — >n c j~^ 3j 

The analysis is based on the study of simulated 
h c events and a background sample of pp — > 37 
candidates (prior to the r\ c mass constraint) from 
a fraction of data taken ouside the Xc.o.g, then 
counting r/ c j candidates in the whole sample. 

Candidates for neutral final states are selected 
by the trigger if there are no charged tracks from 
the interaction point and there are at least 2 en- 
ergy deposits in the Central Calorimeter (CCAL) 
with invariant mass > 2.2 GeV, or if the energy 
detected in CCAL is >80% of the total energy. 

Offline we select events with 3 "on-time" can- 
didate photons in CCAL, defined as clusters with 
energy > 20 MeV occuring within 6 ns of the trig- 
ger. We reject events with on-time clusters in the 
forward calorimeter. Since timing efficiency and 
resolution are worse for low energy photons, clus- 
ters without timing information, or with E < 300 
MeV and occuring within 6 and 15 ns of the trig- 
ger are considered undetermined. If a candidate 
photon paired with any other candidate or unde- 
termined cluster forms a tt° (|m 77 — m^o\ < 35 
MeV) the event is rejected. Events with signals 
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Figure 1. Dalitz plot for candidates pp— >37 events 
in the cos (9* acceptance described in the text 
(small dots) for 3524 < E CM < 3527 MeV 
and pp— *77c7 candidates (large squares). The 
dotted line shows the center of the n c j band. 
pp— >77 c 7candidates with M 2 3 > 1 GeV have rel- 
atively little background. (Preliminary) 



in the two outer hodoscopes not in coincidence 
with the corresponding Cerenkov counter are re- 
jected; those with coincidences are retained as 
events in which a 7 converted after the inner- 
most hodoscope. We impose a likelihood ratio 
test (PW) on photon showers analogous to the 
electron weight (EW) described in Ref. [H], but 
based only on CCAL cluster moments. We re- 
quire PWi x PW 2 > 1 and PW 3 > 1.5, where 
the the photons are ordered by their CM ener- 
gies: Eji > E 7 2 > £73 The efficiency of this cut 
on clean J/ipJ — * e + e~7 events at \ci and \c2, 
of energies comparable to those of h c radiative 
decay, is constant and it is well modeled by the 
simulation. 

A 4C kinematic fit to the hypothesis pp^3j is 




E cu (MeV) 



Figure 2. E835 integrated luminosity vs Ecu in 
the Xc.o.g. region. Dashed histogram is the corre- 
sponding E760 luminosity. Superimposed (dots, 
scale on right) we show the \ 2 probability distri- 
bution vs Ecm for Vcl candidates. Events with 
large probabilities are observed in a narrow Ecm 
range of our high-luminosity scan. Arrow indi- 
cates the probability cut applied. (Preliminary) 



performed and we require a nominal \ 2 proba- 
bility P(37) > 10~ 4 . If there are undetermined 
clusters, we require that P(37) >P(47), the latter 
being the probability for any fit to pp— >4r/. 

37 candidates have background mostly from 
pp— >7r°77, nrf, and Sir events, where the 7r°s 
decay symmetrically or highly asymmetrical; for 
such 7r°'s only one cluster, in approximately the 
direction of the 7r°, is detected and identified as a 
7. We refer to such events as as feed-down events. 
pp annihilations to two and three light mesons are 
strongly forward/backward-peaked and high en- 
ergy photons from feed-down events tend to have 
large |cos#* 2 |. pp^tt°X° reactions are forward- 
backward symmetric, but forward tt° are more 
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3522< E cu <3530 MeV 



3522 3524 3526 3528 3530 
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cand. 


3300 - 3400 






8.51 
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3400 - 3440 






21.83 


3 


3440 - 3500 


0.50 




2.51 




Xcl, 


Xc.o.g reg 


ion and 


Xcl 




3500 - 3520 


5.32 




5.98 




3520.0- 3525.7 


11.24 




13.40 


1 


3525.7- 3526.2 


17.61 


7 


11.56 


6 


3526.2- 3526.7 


5.23 


1 


17.53 




3526.7- 3540.0 


8.12 




8.22 




3540 - 3560 


11.10 


1 


0.90 




3560 - 3675 


33.99 


2 


0.78 




3680 - 3700 


9.00 




12.58 




3700 - 3850 


0.98 




6.32 




3850 - 4400 


1.63 




2.10 





Table 1 

r/c-f candidates in each energy interval and corre- 
sponding integrated luminosities (Preliminary) 



Figure 3. Cross section for Jjp^ry c 7— >777- The 
inset is an expanded view of the Xc.o.g region. 
E835-I data are open circles and E835-II data 
closed triangles. (Preliminary) 



likely to be misidcntificd as single photons, thus 
feed-down events have often low energy 7 s emit- 
ted backwards, with invariant mass of the two 
lowest energy candidates M23 « Mx- 

For reaction (0 the angular distribution is 
nearly uniform in cos 9\ 2 f° r the two photons from 
n c and oc sin 2 #3 for the radiative decay photon. 

Based on signal and background characteris- 
tic we cut at A/23 > lGeV, |cos#* )2 | < 0.5 and 
-0.4 < cos 0*3 < 0.7. 

We apply a 5C kinematical fit to r/ c j and set 
a x 2 probability cut at > 0.01. The overall effi- 
ciency of this selection is 3.2%. 

The final sample has 23 T] c j candidates; 15 of 
them are at Xc.o.g (Fig. 01. Candidates in each 
energy interval and the corresponding integrated 
luninosities are listed in Table The observed 
cross section for pp—^rjc^ is plotted in Fig. [3] 
The background is large near Ecm = 3400 MeV 



but decreases rapidly with energy. We performed 
several checks, in particular: we have analyzed 
the two data sets separately and we find com- 
patible excess at the same mass; If we impose 
in the fit M Va = 2850; 3150 MeV, the cross sec- 
tion agrees very well with smooth background; 
For data outside the Xc.o.g. we rescale all energies 
by 3526.2/ Ecm and verify that events would not 
form an n c . 

We estimate the significance of the excess in 
the rjcj channel with several methods: 

1. Binomial significance with a priori hypoth- 
esis: We calculate the cumulative probability for 
detecting at least N s candidates in an a priori 
signal bin, having observed N to t candidates in the 
Ecm range that extends from the \ci to Xc2 , un- 
der the hyothesis of constant cross section. As 
signal bin we take 3525.6< E C m <3526.4 MeV 
where E760 observed an excess of J/ipir . 

2. Binomial significance with a posteriori hy- 
pothesis with correction for multiple hypothe- 
ses: We observe the excess of events in a 0.5 
MeV bin between 3525.7 < E CM < 3526.2 MeV 
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(MeV) 
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(fb) 
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(nominal) 


(sim) 


(sim)* 


0.5(fixed) 


3525.83^Q - ^g 


4 6 +1 - v 

1.5 


77+30 
' ' -24 


-0.36 


34.19 


0.30-10" 3 


0.98-10" 3 


0.86-10" 3 


1.0 (fixed) 


3525.811^2 


5 5+ 23 

°-°-1.9 


oi +32 
8i -25 


-0.38 


35.66 


1.5-10- 3 


3.16-10" 3 


2.76-10" 3 


No resonance: 




-0.74 


40.72 





Table 2 



Fits to the i] c j cross section in the energy range 3300< Ecm <4400 MeV. P (nominal) is calculated from 
X 2 w —2 AlnC 'P(sim) and 'P(sim)* are the probabilities that AlnC or both AlnC and rp p B n< , 7 exceed 
the experimental values on simulated experiments. (Preliminary) 



or a 1.0 MeV bin between 3525.7 < E C m < 
3526.7 MeV. Since this bin is chosen a posteriori, 
the significance is estimated from the cumulative 
binomial probability, calculated as above, multi- 
plied by a conservative factor 10 (or 5) for the 
number of possible signal bin choices. 

3. Poisson significance: From a linear fit to 
the background cross section o~b over the full en- 
ergy range 3300 < E C m < 4400 MeV, we esti- 
mate crf,(3526.2 MeV) = 0.079 ± 0.016 pb. We 
then estimate the significance from the probabil- 
ity that the expected background fluctuates to 
> 13 events in 3525.7 < E CM < 3.526.2 MeV, 
or > 14 events in 3525.7 < E CM < 3.526.7 MeV, 
multiplied by 10 (or 5). 

4. Significance from likelihood ratio: We sim- 
ulate the outcome of 50,000 experiments under 
the hypothesis of a linear background, whose pa- 
rameters are Gaussian distributed with mean and 
variance taken from the "no resonance" fit to the 
data in Tab. For each experiment we perform 
maximum likelihood fits to the null hypothesis 
(HO) (no resonance) and the alternate hypoth- 
esis (HI) that includes a resonance as described 
below. We then estimate the significance from the 
probability that a likelihood ratio at least as large 
as that observed on data can arise by chance. 

The most conservative estimate of the signif- 
icance is obtained by method 4 which gives V 
between 1 and 3xl0 -3 depending upon the as- 
sumed resonance width. Other methods yield 
8T0~ 5 < V < 3T0~ 3 . In the absence of a narrow 
peaking background, this is strong evidence for a 
resonance near 3526 MeV. 

We perform a Poisson maximum likelihood fit 
to the measured cross section between 3300 and 
4400 MeV as the sum of a linearly varying back- 



ground cross section (crb(E) = cto + b(E(MeV) — 
3526.2)) and a Breit Wigner convolved with a 
Gaussian describing the beam energy distribu- 
tion. 

The parameters determined by the fit are ao, 
b, Mr, and T in B out = T(h c — > pp)B(h c — > 
r] c ^j)B(r\ c — > 77). Data are insufficient to fit for 
r# and we perform fits for fixed values of be- 
tween 0.5 and 1 MeV. The results are given in 
Table for two values of Tr . The background 
parameters are relatively independent of Tr, and 
TinBout only changes by « 20% as Tr is icreased 
from 0.5 to 1.0 MeV. The maximum of the likeli- 
hood seem to favour smaller T^. 

We also perform a fit including in the model the 
oc sin 2 #3 distribution of the El radiative transi- 
tion. The background angular distribution, mea- 
sured on a sample enriched with feed-down events 
at the Xc.o.g. region, is compatible with isotropy 
in —0.4 < cos #3 < 0.7. The fitted parameters do 
not significantly change but InC increases by 2.68 
suggesting that the expected angular distribution 
is a better hypothesis than isotropic decay. 

Dividing T in B out by the value of B(r\ c — > 77) = 
(4.3±1.5)10" 4 g|, we derive 10.6±3.7±3.4(6r) < 
Tp P B Vcl < 12.8±4.8±4.5(6r)eV, consistent with 
B(pp) w 1 — 3xl0~ 5 in the expected range [TU| 
forr(/i c ->7? c 7)« 600keV. 

3.2. Search for pp — > h c — > J/tpir — * e + e _ 77 

The trigger for e + e~ X final states requires 
at least two charged tracks from the interac- 
tion point associated to a signal in the threshold 
Cerenkov counter and at least 2 energy deposits 
in the CCAL with invariant mass > 2.2 GcV. 

Events must have two electrons identified by a 
likelihood ratio test (EW X EW 2 > 1.5 0) based 
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Figure 4. Cross section for pp— ► J / ipir e~ 'yy . 
E835-I data are open circles and E835-II data 
closed circles (Preliminary) 



on dE/dX in the scintillators, number of ph.c. 
in Cercnkov and shower lateral shape in the 
CCAL. We require that the e + e~ invariant mass 
is > 2800 MeV and limit the e + (e~) acceptance 
to 15° < 8 < 60°. Acceptance for photons is 
11° < 6 < 70°. We allow additional on-time 
CCAL clusters only if compatible with photons 
radiated by the e ± (E < 100 MeV and 6> 7e < 10°). 
Finally the x 2 probability for the 6C fit to J/ipir 
must exceed 0.01. 

We exclude the \a and %c2 data, since ra- 
diative decays to J/i/jj constitute a background 
to J/ipn°, and consider only data for 3.52 GeV 
< < 3.54 GeV. The observed cross section for 
pp^J/if)ir° is plotted in Fig.0] Data are compati- 
ble with a flat cross section between 3.52 and 3.54 
GeV. Within our acceptance (smaller than that 
of E760) we see no excess of events that would 
correspond to a narrow (< lMeV) resonance. 

4. CONCLUSIONS 

We have measured the cross section for pp — ► 
Vcl ~ * 37 near the center of gravity of the char- 
monium 3 Pj states and observe a narrow struc- 



ture (r < 1 MeV) centered at 3528.8 ± 0.2 ± 

0. 2 MeV. The statistical significance is V ~ 0.001. 
The value of V Pp B Vcl ~ 12.0±4.5±4.3eVis com- 
patible with the expected B(pp) and El radiative 
width of the h c . The signal is seen with com- 
parable yield and at the same Mr in both runs. 
Our former experiment E760 had no sensitivity 
to this channel: with a factor 5 less luminosity in 
this region, it would have observed between 2 and 
3 events. This observation is evidence for the h c 
at its expected location near the Xc.o.g- We have 
measured the cross section for pp —> J/ipir and 
observe no significant excess in the Xcj center-of- 
gravity region. 
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